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I, Introduction 



A. Overview 

The purpose of this report is to discuss solution methods for several 
problems arising in inventory provisioning. Two related types of problems are 
considered: 

1) optimization of a performance measure subject to a budget 
constrai nt 

2) minimization of cost subject to a constraint on performance. 

The performance measures considered are Supply Material Availability (SMA), 
Mean Supply Response Time (MSRT), and "Pseudo-Availability" (PA). The basic 
approach to these problems is dynamic programming. 

First the problems are formulated and solved using one recursive 
technique. Then a more efficient recursion is presented and discussed in the 
context of maximizing MSRT subject to a budget constraint. 

These same budget constrained problems are formulated, discussed and 
solved using a marginal analysis approach in [2] and [3], That method, 
although fast, does not guarantee that optimal solutions are obtained. This 
report provides a method of obtaining optimal solutions and thus provides a 
means of evaluating heuristic methods. In addition, this report shows how the 
performance constrained budget minimization problems can be solved directly. 

By contrast, reference [2] addresses this problem using generalized Lagrange 
Multipliers or by solving the budget constrained performance problem 
repeatedly for various budget levels. 

The second recursion presented here is considerably faster than the first 
and it also guarantees an optimal solution in some cases. It is competitive 
in speed with the marginal analysis method for small and medium-sized 
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problems, but it is inefficient for problems with a large number of items or 
large budget values. Its virtue lies in its ability to get exact solutions to 
medium-sized problems. 



B. Problem Formulations , budget constrained 

This report considers three specific budget constrained optimization 
problems arising in inventory provisioning. These problems are formulated and 
discussed in detail elsewhere [2] and will be stated here without extensive 
explanation. 

Following [2] we let 



n 


= the total number of items considered for provisioning 


'i 


= the unit cost of item i 


(- 

n . 


= the essentiality code for item type i 




= the demand rate for item i 


T. 

1 


= the procurement leadtime for item i 


s. 

1 


= the number of items of type i provided (the decision 
vari ables ) 


Z,(S,) 


= the performance measure for item i when S. units are 
stocked 


0,(S.) 


= - Z,(S,) 


P,(x,) 


= probability that demand for item i is x. during the 
provisioning interval 


P,(x,) 


= cumulative probability of x^ or fewer demands during 
a provisioning interval 


MTTR. 

1 


= mean time to repair or replace item i 


MTBF . 
1 


= mean time between failures = 1/X^. 


MSRT.{S. ) 
1 ' 


= mean supply response time when S^. units are stocked. 
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where 



where 



where 



The three budget constrained problems considered are: 

al) maximize Supply Material Availability (SMA), defined as 



SMA(S^,... ,S^: 



n 

= I 

i=l 



E. 

1 



X. 

1 



T. 

1 



)/ I 

^ ’ i=l 



E. 

1 



X. 

1 



T. 

1 



(1) 



zr^(s.) = (1 - p.(s.)) + (S. - X. T.)(l - P.(S.))/X. T. 



bl) minimize Mean Supply Response Time (MSRT), defined as 



MSRT(S,,...,S ) = y E. X. T. zf^^S.)/ I E. X. T. 
W’ " n' 1 1 1 1 ^ ^ 1 1 



zl^^S. ) = (1 

1 ' V 



P.(S.))(T./2 - S./X. + 
1 ^ 1 . ^ 1 11 



2X^ T. 
1 1 



+ p.(S.)(X. T. - S.)/2X. 
*^111 1 1 ^ 



cl) maximize Pseudo-Availability (PA), defined as 



PA(S,,...,S ) = II zf^S 
i n 1 



zP^s 



= MTBF./(MTBF. 



+ MTTR. 
1 



+ z!^^s. )) 

1 1 
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In each of the budget constrained problems above the "objecti ve is to 
be optimized by selection of subject to the constraints 



where B is the specified budget level. 

C. Problem Formulations, Performance Constrained 

In addition to the three problems just stated, we consider three 
related problems in v/hich the cost is to be minimized subject to a constraint 
on performance. 



n 



I i. ^ 1 ^ integer 



n 



a2 ) min T C . S . 

i=l ^ ^ 



s.t. SMA(Sp 



S^) ^ SMA 



S . >_ 0 i nteger. 



n 



b2) min I C. S. 

i=l ^ ^ 



s.t. MSRT(Sp 



• • • J 



S) < MSRT 
n — 



S . ^ 0 i nteger. 



n 



c2) min I C . S . 

i=l ’ ^ 



s. t. PA(Sp... ,S^) >_ PA 



^ 0 integer. 
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II. Solution Method and Examples 

A. Dynamic Programming Approach 

The computer program used to solve these problems is DP4, a general 
purpose program for performing dynamic programming tabular computations. Here 
we will describe the general nature of that program and those elements 
required to tailor it for use in the problems considered in this report. 

The DP4 program deals with a problem consisting of n related stages 
each of which is characterized as shown in figure 1. 



X. 



U . 

i: 



T 






t,(x,,d.) 



r. (X. ,d. ) 

V 1 ’ V 

Figure 1. A single-stage decision problem. 



In figure 1 



X- is the "state" variable 

d^. is the decision variable 

r^ is the stage return function 

t. is the stage transformation function. 

In the overall problem consisting of n stages the output state variable from 
stage i , namely X^._^ , is the input to stage i-1 . Thus, the n stage 
problem can be pictured as in figure 2. 



n-1 






d. 

1 






Figure 2. n-stage decision problem. 
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The state variable X^. can easily be understood in the context of 
the budget constrained problems as the amount of resource (money) remaining to 
be allocated to stages i, i-l,...,l . The stages, of course, correspond to 
the items in the inventory problems. 

At each stage i a decision d^ must be made. The decision has two 

effects. First, it yields a return r^. , the performance measure associated 

with the current item. Second, it yields a value of X^. ^ which serves as 

the input to the remainder of the decision process. The decision d^ must be 

made with consideration both for the immediate return r^ and the future 
state X^_^ . The overall problem is to make the series of decisions 
d^,...,dj^ to optimize some function of the infdividual stage returns. 

In the problem (al), where S^. is the decision variable, we can let 
the return functions be 



r,(X., S,) . E, X, 1. Z<*>(sp/,I^E, X, T, 1 = 1 n 



and the stage transformation functions be 






The overall return function is the sum of the individual return 
functions. Namely, 



SMA(S ,...,$ ) = I r.(X , S.) . 

1 n 1 1 1 

The object remains to select S^,...,S^ to optimize this return. 
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We let = the optimal total return from stages i, 

given that we enter stage i with state variable X.. 
Then we can write the recursive equations for this optimization as 

*’((X() ' ™x{r,(X,, s.) t 

s.t. X. , = X. - C. S. 

1-1 1 11 

and 0 < S. < X./C. 

- 1 — 11 

and = integer 

for i = 2,. . . ,n . 

The equation at stage 1 is 

f^(X^) = max n^(Xp S^) 

s. t. 0 _< _< 

and integer . 

The program DP4 performs this optimization provided the user supplies 
the following subroutines and data. 

Required subroutines 

1. STGRET - this subroutine defines the function '^■j ) 

2. TRANFM - defines the stage transformation function t.(X., d^ ) 

3. OLIMIT - defines the range of decision values d^ which can be 

considered for the particular value of X^ under 
consi derati on 

4. STORE - allows the input of constants to be used in the other 

subrouti nes. 
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Required Data 

1. n - the number of stages 

2. For each stage i 

XLOW - the lowest value of to consider 

XHIGH - the highest value of X^ to consider 

OELX - the increment for X^ 

XMODE - tells whether to maximize or minimize 

XSTAGE - tells how this stage return relates to lower numbered stage 
returns (sum, product). 

The methodology is essentially the same for the performance con- 
strained problems. There the return functions '^■j ) " • The 

state variable X^ is interpreted as the portion of the performance measure 

to be attributed to stages 1 i . The stage transformation functions in 

problem (a2) and (b2) are 

*(-l ■ • 

In problem (c2) the stage transformation is 

Tne program DP4 and the subroutines are shown in Appendix A. The 
subroutines are written to solve any of the problems (al), (bl), (cl) or (a2), 
(b2), (c2). Thus they involve complications not needed for solving just one 
of these problems. 

B. Examples 

Several example problems were solved to illustrate the approach discussed 
here. All of the problems involved n = 10 items and all used the data shown 
in table 1. 
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Data 



n 


X 


Time 


Cost 


MTTR 


E 


1 


b.O 


1.0 


1.0 


.0137 


1.0 


2 


2.0 


1.0 


2.0 


.0274 


1.0 


3 


3.0 


1.0 


b.O 


.0137 


1.0 


4 


b.O 


1.0 


10.0 


.0822 


1.0 


5 


10.0 


1.0 


20.0 


.0274 


1.0 


6 


2b. 0 


1.0 


b.O 


.0027 


1.0 


7 


1.0 


1.0 


1.0 


.0054 


1.0 


8 


1.0 


1.0 


lOO-.O 


.0411 


3.0 


9 


O.b 


1.0 


bO.O 


.0082 


1.0 


10 


2.0 


1.0 


10.0 


.1370 


3.0 



Table 1: Data for Examples 

1. The budget constrained problems. 

Tables 2, 3, and 4 summarize the solutions for the example problems (al), 
(bl), and (cl). These are the budget constrained problems. 
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B = Budget = 
max SMA = 


300 

.750654 


295 

.746466 


290 

.741160 


285 

.736563 


280 

.732964 


decision = 


4 


4 


4 


4 


4 


^2 - 


4 


4 


4 


4 


4 




4 


4 


4 


4 


4 


^4 - 


5 


5 


5 


4 


5 




2 


2 


2 


2 


1 




29 


23 


27 


28 


29 




3 


3 


3 


3 


3 




0 


0 


0 


0 


0 


Sg = 


0 


0 


0 


0 


0 


^10 " 


3 


3 


3 


3 


3 



Table 2: Solutions to example problem (al) 
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B = 

min MSRT = 


300 

.0838231 


295 

.0859625 


290 

.0878466 


285 

.0900560 


280 

.0927486 


decision 


2 


2 


2 


2 


2 


^2 


3 


3 


3 


3 


3 


^3 


3 


4 


3 


3 


3 


S 4 


4 


4 


4 


4 


4 




6 


4 


4 


4 


4 


^6 


21 


23 


23 


22 


21 


^7 


2 


2 


2 


2 


2 


Sb 


0 


0 


0 


0 


1 

0 


S 9 


0 


0 


0 


0 


0 


^10 


3 


3 


3 


3 


3 



Table 3. Solutions to example problem (bl) 
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8 = 
max PA = 


300 

.0403727 


295 

.0387421 


290 

.037005 


285 

.0355446 


280 

.0341911 


decision S^ 


3 


3 


3 


3 


3 


^2 


4 


4 


4 


4 


4 


^3 


5 


4 


5 


5 


5 


^4 


t) 


5 


5 


5 


5 


^5 


3 


3 


2 


2 


2 


^6 


26 


26 


26 


27 


26 


^7 


4 


4 


4 


4 


4 


^8 


0 


0 


0 


0 


0 


^9 


0 


0 


0 


0 


0 


^10 j 


2 


2 


3 


2 


2 



Table 4. Solutions to example problem (cl) 

2. The performance constrained problems. 

The related performance constrained problems were also solved for illus- 
tration. For example, the problem (a2) was solved using the data from Table 1 
with the restriction that SMA >_ . 732964 . The solution to that problem is the 
same as the solution shown in the last column of Table 2 since the value 
.732964 is the (largest) value of SMA corresponding to a budget of 28U. 
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III. Modifications to Basic Method 



A. An Alternative Recursion 

An alternative and more efficient approach is available for the budget 
constrained problems. 

The approach is based on a different recursion from that used in 
Section 2 and is very similar to an approach used to solve the "cargo loading 
problem". See for example Dreyfus [1]. 

The cargo loading problem, stated as a maximization, is: 

N 

max y V .d . 

j=i J J 

s.t. UjCljlB 

d j _> U i nteger . 

Although many methods are available for solving this problem, the one of 
interest to us is based on the following recursion 

f(b)= max {v.+f(b-c.)} 
je{l,...,N} ^ ^ 

where f(b) is the optimal total return that can be obtained when a budget 
of b is avai lable. 

To illustrate this method consider the data in Table b. 



i 


1 


2 


3 = N 




o < 


1 


4 


6 


B = 10 


1 


3 


4 





Table b. Data for example using alternative recursion 

The computation proceeds with increasing values of b until b = B is 
reached. The process can be viewed as shown in figure 3 where a template 
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representing the available items is placed over the budget value of current 
interest. The template points back to previously determined optimal solution. 
Each of these previous solution is considered for updating by including one 
more item of the type inducted by the template. The updated solutions are 

compared and the best selected as the solution for the current value of b . 

The illustration shows the template at the budget value of 7. The optimal 
solutions for b = 0 , 1 ,..., 6 have already been computed. The comparison at 
B = 7 is among the solution at 6 with an additional item 1 for a total return 
of 9, the solution at 4 with an additional item of type 2 for a total return 
of 10, and the solution at 3 with an additional item of type 3 for a return 

of 10. Either of the last two is chosen and recorded as an optimal solution 

at 7. 




b = 0 



2 3 4 5 6 



8 9 10 



d.(b) = 0 

dp(b) = 0 

d3(b) = 0 



f(B) = 0 



1 2 0 0 1 
0 0 10 0 
0 0 0 1 1 



2 * 0 . 

0 1 

1 1 



12 4 6 



8 10 



*an alternative optimal 
solution at b = 6 is 
d^ = 0 d2 = 2 d3 = 0 



Figure 3: Illustration of solution method using the alternative recursion. 



A very simple modification of this procedure can be used to solve the 
problems of the type discussed in this report. We consider 

n 

max I r.(d.) 

j = l J J 



,t. I c . • d. < 
j = l ^ ^ 



d. > 0 integer 
J ~ 
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where the return functions r.(d.) are concave. The procedure below is not 

\J 

guaranteed to give the optimal solution for all r.(d.) but is guaranteed 
if the r.(d.) are points on a concave function. 

\J \J 

In this case we will represent the return functions 
sum of the marginal values of additional items of type j . 

d 

r .(d) = + I m. (i ) 

J i=l ^ 



Thus r.(2) = m.(U) + m.(l) + m.(2) . These marginal values form a 

\J \J U \J 

sequence with the properties that 



and 



m . ( i ) > 0 
J 1 

m^ (i ) > m^ (k ) i < k . 



The same algorithm as before was applied with the modification that the value 
term v- , which was formerly constant, is replaced by m-(i) for the 

\J \J 

appropriate value of i . The program which implements tnis algorithm was 
called RECUR. 

It should be noted that the discussion above treats the constraint as 
an inequality, but the function f(b) in this section is computed for the 
constraint 



N 

V 

L 



j = l 



C .X . 



J J 



= b . 



For this reason we may have in a maximization problem 

f(b^) > f(b2) 
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although < b 2 . That is, if we require the equality to be met exactly, 
it is not necessarily true that a larger budget is better. The program prints 
the values of f(b) for several values of b so the optimal value of b _< B 
can be found visually. 

B. Minimum Orders 

1) RECMOD 

A modification was made to the program RECUR to permit the user 
to specify minimum packaging quantities of each item. That is, item i is 
assumed to be packaged with q^ items per package. The provisioning can 
select only whole packages of each item. This modification resulted in the 
program RECMOD given in the Appendix. 

2) Example 

To illustrate the RECMOD program, consider the example problem 
(a2) solved previously. The optimal solution for budget B = 300 is repeated 
in Table 2 for the case in which all q^ = 1 . 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Object i ve 


^i 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




^i 


2 


3 


3 


4 


6 


21 


2 


0 


0 


3 


.0838231 


^i 


2 


1 


1 


5 


1 


1 


1 


3 


1 


1 




^i 


2 


3 


3 


5 


4 


23 


2 


U 


0 


3 


.0847410 



Table b. Solution to Example using RECMOD. 



The optimal solution is also shown for a modified problem in which not all q^- 
are equal to 1 . 
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C. Discussion 



The modified recursion just discussed has been implemented for the prob- 
lem (a2) which is to minimize MSRT subject to a budget constraint. The method 
can also be applied to the other budget constrained problems, but this has not 
yet been done. All that is required is to modify the program to compute SMA 
or PA instead of MSRT and to maximize instead of minimize. 

There is a difficulty in extending this method to problems in which the 
item costs are arbitrary values. The method is very effective when the costs 
are all integer and can be scaled so that the smallest cost is 1 . If 
arbitrary costs are allowed, the algorithm can become ineffective for all 
except small values of B . Consider for example the costs of Sl.UU, $1.21, 
$1.27 for three items. Let the budget be $2b.U0. The problem could be solved 
by scaling the costs to be 100, 120, and 125 and the budget to be 25U0, but 
then too many values of b must be considered when many of them are not pos- 
sible. Alternatively the algorithm can step to the "next possible value" and 
will consider the folowing sequence of values 

b = 100, 121, 127, 200, 221, 227, 242, 248, 254, 300, 321, 327,... . 

As the process continues, depending on the relative values of the costs, the 
sequence becomes more dense and may eventually include all possible values of 
b . This is ineffective and cumbersome for large values of B . 

It is also not possible to apply the modified recursion for the perform- 
ance constrained problems. On the other hand, the budget constrained problems 
are solved very rapidly and the relationship between cost and performance can 
easily be determined from solving the budget constrained problem using the 
first recursion. In fact, if MSRT is minimized for a budget of B , the 
solution is actually obtained for all values of b up to and including B . 
Those results reveal the relationship between performance and budget. 
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APPENDIX-PROGRAM LISTINGS 



FILE: DPFILE EXES A1 NAVAL POST3RADJATE SCHOOL 



FILEDEF 07 DISK 
FILEDEF 08 DISK 
FILEDEF 06 TERM 
FILEDEF 02 DISK 



(PERM 

RUSSIO DATA (PERM 
(LRECL 133 RECFM F8 PERM 
( PERM 



FILE: FILE 



FT07F001 A1 



NAVAL POSTGRADUATE SCHOOL 



10 


NEW 






1 


10 


0. 


3 00, 




300. 




300, 


10 


CLO 








2 95. 




2 9 5. 


10 


CLD 








290. 




2 90, 


10 


OLD 








285. 




2 35, 


10 


OLD 








280. 




2 80, 


10 


CLD 








275. 




275, 



MIN SUM 2 



FILE; RUSSIO DATA A1 NAVAL POSTGRADUATE SCHOOL 



2 0 

1. 


1. 


1. 


1 . 


.0137 


2* 


2 . 


1. 


1 . 


.0274 


3. 


5# 


1. 


1 . 


. 0137 


5 . 


1 3 • 


1. 


1 . 


.0822 


10. 


20. 


1. 


1. 


. 0274 


25. 


5 . 


1. 


1. 


. 0027 


1 • 


1 • 


1. 


1 . 


.0054 


1. 


100. 


3. 


1. 


. 0411 


.5 


5 3. 


1. 


1 . 


. 0082 


2 . 


10. 


3. 


1 . 


.1370 


10 300(08, 


CON) 3C3ST 


ISMA 2MSRT 


3AVA I L - 


DATA ACEBD = LAM,COST 


2 1 


1 5 


1 1 


1 3 


1 1 


FOR RECUR THE 
CURRENT L INE. 


LAST TWO LINE MUST 
IT IS FOR OPA 


8E PLACED 


FIRST INSTEAD OF TH 



T,M? 
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oooooonoooooonooooooooooooooooooooooooonnoooooooonooooonoooooooooooooooo 



FILE: DP4 



FORTRAN A1 NAVAL P0ST3RA0JATE SCHOOL 



DATA DESCRIPTION FOR CARD NJMBER 1 



COLJMNS 



JUSTI FY 



VARI BLE 
NAME 



MEAN in: 



1-5 



6-10 



C3L 5 



C3L 10 



MMAX 

MTAPE 



THE NUMBER OF STASES IN THE PROBLEM 

L3SICAL TAPE NO OF THE MAIN TAPE 
IF THIS IS LEFT BLANK, THE COMPUTER 
WILL USE STANDARD SCRATCH TAPE 3 



1 1 - 15 



COL 15 



SOLVE 



= OLD IF THE OPTIMAL DECISION 

FUNCTIONS HAVE ALREADY BEEN 
CALCULATED AND ARE LOADED 
ON LOGICAL TAPE MTAPE 



= NEW IF THE OPTIMAL DECISION 
FUNCTIONS MUST BE 
CALCULATED AND STORED ON 
LOGICAL TAPE MTAPE BEFORE 
SOLVING THE PROBLEM 



DATA DESCRIPTION FOR THE NEXT GROUP OF CARDS 
IF SOLVE = NEW ON CARD NUMBER 1, THEN EACH OF 
MUST BE DECRIBED BY THE CARDS DISCUSSED BELOW. 
PER STAGE OR YOU CAN MAKE 1 CARD DESCRIBE MANY 
THEY ARE SIMILAR. THE STAGES MUST BE DESCRIBED 
STARTING WITH STAGE 1 (I.E. STAGE 1, THEN 



THE iNMAX STAGES 
YOU CAN USE 1 CARD 
ADJACENT STAGES IF 
IN NUMERICAL ORDER 
THEN NMAXJ 



OIM’^IXINO olM (jC i. I X . C . O I M j Z i. , I nC I't ..., inCIN 1,11 MAI 

OMIT THIS PACK OF CARDS WHICH DESCRIBE THE STAGES IF SOLVE = OLD 



COLUMNS 


• 

• 


justify I 

• 


VARI BLE 
NAME 


1 MEANING 


1 - 


5 


* • 
• 
• 


COL 


• 

5 

• 


NSTAGE 


I LOWEST NUMBERED STAGE FOR WHICH 
. THI S CARO APPLI ES 


6 - 


10 


• 

• 

• 

• 


COL 


10 I 

• 

• 

• 


NDITTO 


I HIGHEST NUMBERED STAGE FOR WHICH 
. THIS CARD APPLIES. IF THESE COLUMNS 
. ARE LEFT BLANK, THEN NDITTO WILL BE 
. TAKEN AS EQUAL TO NSTAGE 


11 - 


20 


• 

• 


ANY 


• 

• 


XLOW 


I LOWEST VALUE OF XN FOR STAGE NSTAGE 


21 - 


30 


• ■ 
• 


ANY 


• 


XHIGH 


1 HIGHEST VALUE OF XN FOR STAGE NSTAGE 


31 - 


40 


• 

• 


ANY 


• 


DELX 


I INCREMENT IN XN FOR STAGE NSTAGE 


41 - 


46 


• 

• 

• 


COL 


46 1 

• 


XMODE 


I = MIN IF STAGE NSTAGE IS TO BE 

. MINIMIZED 






• 

• 

• 




• 

• 

• 




I = MAX IF STAGE NSTAGE IS TO BE 

. MAXIMIZED 


47 - 


52 


• 

• 

• 

• 

• 


COL 


52 I 

• 

• 

• 


XSTAGE 


r = SUM IF THE COMPOSITION 

. OPERATOR BETWEEN STAGES 

. NSTAGE AND NSTAGE-1 IS 

. ADDITION 






• 

• 

• 

• 

• 




• 

• 

• 

• 

• 




1 = MULT IF THE COMPOSITION 
. OPERATOR BETWEEN STAGES 

. NSTAGE AND NSTAGE-1 IS 

. MULTIPLICATION 



=M INMAX IF THE COMPOSITION 

OPERATOR BETWEEN STAGES 
NSTAGE AND NSTAGE-1 IS TO 
MINIMIZE THE MAXIMUM 
INDIVIDUAL STAGE RETURN 
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53 



=^AXMIN 



55 



COL 55 



MODES 



= 1 



IF THE COMPOSITION 
OPERATOR BETWEEN STAGES 
NSTA3E AMO NSTAGE-1 IS TO 
MAXIMIZE THE MINIMUM 
INDIVIDUAL STAGE RETURN 



THIS VARIABLE IS NOT USED 
WHEN NSTAGE = 1, THEREFORE 
COLUMNS '■*! TO 52 ARE 
IGNORED WHEN NSTAGE=1 AND 
MAY BE LEFT BLANK 
HOWEVER IF NDITTO IS 
GREATER THAN 1, THEN XSTAGE 
MUST BE SPECIFIED IN 
ACCORDANCE WITH THE 
COMPOSIT3N OPERATOR FOR THE 
REST OF STAGES DESCRIBED ON 
THIS CARO 



IF THE FIBONACCI SEARCH IS 
TO BE USED 



IF THE FIBONACCI SEARCH IS 
NOT TO BE USED 



DATA DESCRIPTION OF LAST CARO FOR PROBLEM 

THE OPTIMUM DECISIONS ARE TO BE PRINTED OUT FOR A PROCESS WITH 
NMAX STAGES, FURTHER, FN(XN» IS TO BE OPTIMIZED WITH RESPECT TOXN 
WHEN THE PROGRAM COMPLETES THE REQUESTED PRINTOUT IT WILL RETURN 
TO THE BEGINNING AND START A NEW PROBLEM IF THERE IS DATA, IF NO 
DATA , IT WILL EXI T 



COLUMNS 


“ • 
• 
• 


UUSTIFY 1 

• 


VARI 3LE 
NAME 


• 

• 


MEANING 








1-15 


• • 
• 
• 


ANY 1 

• 


XNl 


* • “ 

• 

• 


AT STAGE NMAX, OPTIMIZE 
XN NOT LESS THAN XNl 


FN( XN) 


FCR 


16 - 30 


• 

• 

• 

• 


ANY 1 

• 


XN2 


• 

• 

• 


AT STAGE NMAX, 
XN NOT GREATER 


OPTIMIZE 
THAN XN2 


FN(XN) 


FCR 



DESCRIPTION OF THE SUBROUTINE TRANFM 

THE MAIN PROGRAM TRANSMITS XN,ON, AND NUMBER (THE SEQUENCE NUMBER 
OF THE PROBLEM! TO THE SUBROUTINE AND THE SUBROUTINE CALCULATES 
YN (THE OUTPUT OF THE STAGE I.E. YN=XN-1) 



DESCRIPTION OF THE SUBROUTINE STGRET 

THE MAIN PROGRAM TRANSMITS XN,DN,YN, AND NUMBER TO THE SUBROUTINE, 
AND THE SUBROUTINE CALCULATES RN (THE STAGE RETURN! 



DESCRIPTION OF SUBROUTINE STORE 

THE MAIN PROGRAM CALLS THE SUBROUTINE RIGHT AFTER READING THE 
FIRST DATA CARD AND DEFINING NUMBER. THE PURPOSE OF THIS 
SUBROUTINE IS TO ALLOW THE STORING OF CONSTANTS IN COMMON STORAGE 
FOR THE POSSIBLE USE OF SUBROUTINES STGRET AND TRANFM 

DATA WHICH WILL BE READ FRCM DATA CAROS AT 3BUECT TIME BY 

THE SUBROUTINE STORE MUST BE INSERTED BETWEEN THE FIRST DATA CARO 

AND THE GROUP OF CAROS WHICH DESCRIBE THE INDIVIDUAL STAGES 
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DESCRIPTION DF SUBROUTINE DLIMIT 

THE MAIN PR03RAM TRANSMITS N.XN, AND NUMBER TO THE SUBROUTINE, 

AND THE SUBROUTINE CALCULATES DL3W, DHIGH, AND OELD 
WHERE DLOW * LOWEST VALUE OF DN FOR THAT PARTICULAR VALUE OF N AND 
XN 

HIGHEST VALUE OF ON FOR THAT PARTICULAR VALUE OF N 
AND XN 

INCREMENT IN DN FOR THAT PARTICULAR VALUE OF N ANC XN 
DO NOT CHANGE, THEY CAN BE SET IN 
MENTIONED IN SUBROUTINE DLIMIT 



DFIGH = 

OELD = 

IF ANY OF THESE VARIBLES 
SUBROUTINE STORE AND NOT 



IF ANY SUBROUTINE SETS NUFF TO -1 THEN THE PROGRAM IGNORES THE 
VALUE OF DN BEING PROCESSED AT THAT MOMENT, AND IGNORES ALL LARGER 
VALUES OF DN FOR THAT PARTICULAR VALUE OF XN AT THAT PARTICULAR 
STAGE — THIS FEATURE IS DISABLED NHEN USING THE FIBONACCI SEARCH 



NUMBER 



0 .0,- 1.0,0.U,DN3E ST,BEST ) 



NMAX,MTAPE, SOLVE 



THIS PROGRAM USES TAPES 3 AND 4 FOR SCRATCH 
NSS2 ZERO FOR ENTIRE PROBLEM SO-VING AND OUTPUT 
NSS2 NON-ZERO FOR OUTPUT FROM PREVIOUS PROBLEM 
DIMENSION FN( lOOOli ,FNM1 ( 1 0001 », 0 NOFX N ( 1 0001 ) 

REAL 8 HOLM IN/' MIN'/, HOLQLO/' OLD'/ 

REAL * 8 XMOOE, XSTA3E, SOLVE 
REAL * 8 TAB(4» 

COMMON /CON/TAB 

COMMON XN,N, ON, YN ,RN, SLOW ,XHIGH, OE LX, DLOW,DHI GH, 0 ELO, 

COMMON YLOW, YHIGH ,DELY,NUFF,FN,UTOP,NNFlAG,NMAX,XLAM 
COMMON A ,B,C, D, E , SUMELT, KODEl, K00E2 
EQUIVALENCE ( FN( 1 ),FNM1( 1 » ,DNOFXN ( 1)) 

WRITE(6, 10001 HOLMIN, HOLOLD, (TAB(I), 1 = 1,4) 

1000 FORMATdX, 6A6) 

MAX= 10001 
INTAPE=7 
N0UTPE=6 
CALL SEARCHO.O 
NPAGE=0 
NUMB ER=0 

240 READ ( INTAPE, 100 ,END=2 000 ) 

100 F0PMAT(2I 5,A5) 

MTAP = MTAPE 
NDITT0=0 
NCALC=4 

IF(MTAPE-4) 405,501,405 

501 NCALC=3 
GO TO 172 

405 IF(MTAPE) 322,322 ,172 

322 MTAPE=3 

172 NTAPE=MTAPE 
REWIND MTAPE 
NUMBER=NUMBER+1 

; CALL STORE WAS MOVED FROM HERE 

IF(SOLVE-HOLOLD) 301,302,301 

302 NSS2=1 
NPAG E = NPAGE+1 
WRITE! NOUTPE , 152 ) NPAGE 

152 FORM AT (1H1100X4HPAGEI 4) 

WRITE! NOUTPE, 303) MTAPE 

303 FORMAT !58H THE TABLES OF 
1APEI5,35H WILL BE USED TO 

READ (MTAPE) NMAX 
REWIND MTAPE 
GO TO 304 

301 NSS2=0 

CALL STORE 
NPAGE=NPAGE+1 
WRITE(NOUTPE, 152) NPAGE 
WRITE! NOUTPE, 305 ) MTAPE 

305 F0RMAT(39H THE COMPUTER IS TO CALCULATE TABLES OF FN(XN) AND ON(XN 



FN(XN) AND 
SOLVE THIS 



DN!XN) WHI 
PROBLEM) 



.H ARE ON LOGICAL T 
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500 

304 

306 

314 

109 

2006 

121 



241 

124 



2007 

31 

30 



600 

601 

125 

127 

245 

246 

247 



), AND STORE THEM ON LOGICAL TAPEI5) 

REWIND NCALC 
WRITE(PTAPE» NMAX 
WRITE(NOUTPEi 306) NMAX 

F0RMAT(16H THE PROBLEM HASI5,7H STAGES) 

IF(NSS2) 197,109,197 
LINES = 6 

WRITE! NOUTPE, 2006 J NMAX, MTAP, SOLVE 
FORMAT!/// • DATA CARD •, 217, A7) 

DO 123 N=1,NMAX 
NM1=N-1 

IF!NDITTO-N» 241,132,132 

READ! INTAPE, 124) NSTA GE, ND ITTO ,X L OW, X HIGH , OE L X, X M ODE , XSTA GE, MODE S 
F0RMAT!2I5,3E10.6,2A6 ,13 » 

NDITTO, XLOW, XHIGH, DELX, XMODE, 



001 



OF ORDER) 



1001 

126 

231 

2 32 
230 

234 

520 

129 

521 

130 
502 
316 
110 

1 

313 

510 

321 

511 

3 20 
11 
16 

1 

12 

17 



WRIT£!NOUTP£, 2007 ) NSTAGE, 

XSTAGE, MODES 

FORMAT!/// ' DATA CARD ', 
IF!LINES-41) 30,30,31 
NPAGE=NPAGE+1 
WRITE!NOUTPE, 152) NPAGE 
LINES=1 

LINES = LINES + 9 
ITOP=! XHIGH-XLOW ) /DELX+1 
IFIMOOES-1) 600,601 ,600 
MODE S=-l 

IF!NSTAGE-N) 125,245, 125 
WRITE!NOUTPE,127) 

F0RMAT!23H1DATA CARD OUT 
GO TO 2003 

I F!NSTAG£-N0ITT0 ) 246 ,126,126 
WRITEiNOUTPE, 247 ) N, N STA G E ,NDI H 3 
FORM AT!///46H THE DESCRIPTION WHICH 
, APPLIES TO STAGE, 15, IIH THRU 
L INES=LINES+4 
XXMODE = XMOOE - HOLM IN 
WRITE! 6,1001) XXMODE 
FORM AT! IX, A6 ) 

I F!XMOOE-HOLMIN ) 230,231 

XMOOE=1.0 

WRITE!NOUTPE,232 ) N 
F0RMAT!//6H STAGEI5,18H 
GO TO 520 
XMQDE=-1.0 
WRITE!NOUTPE, 234) N 
F0RMAT!//6H STAGE 15, 18H 
DO 129 1=1,4 

129,13 0, 



215, 2X, 3Fi0.6, 2Ab, 13) 



APPEARS BELOW FOR 
STAGE,I 5,10H INCLUSIVE) 



STAGE, 15, 17H 



230 



I S A MINI MIZATION) 



IS A MAXIMIZATION) 



129 



TREATED AS DISCRETE, 
STEPS OF E14.5) 



GOING F 



13 

18 



I F!XSTAGE-TAB! I ) 

CONT INUE 
NNFL AG=5 
GO TO 502 
NNFLAG=I 

IF!ITOP-MAX) 316,316,112 
WRIT E!NOUTPE, 110) XLO W ,X H I GH , DEL X 
F0RMAT!63H THE STATE VARIBLE IS BEING 
ROM XN=E14.5,7H TO XN = E14.5,13H IN 
IF!MODES) 318,318,321 
WRIT E!NOUTPE, 510 ) 

FQRMAT!41H TOTAL ENUMERATION IS USED FOR THIS STAGE) 

GO TO 320 
WRITE!NOUTPE, 511 ) 

F0RMAT!40H FIBONACCI SEARCH IS USED FOR THIS STAGE) 

GOTO !11,12, 13, 14,132),NNFLAG 
WRITE1N0UTPE,16 ) NH1,N 

F0RMAT!40H THE COMPOSITION OPERATOR BETWEEN STAGESI6,4H ANDI6,12H 
IS ADDITION) 

GO TO 132 

WRITE!NOUTPE, 17) NM1,N 

F0RMAT!40H THE COMPOSITION OPERATOR BETWEEN STAGESI6,4H ANDI6,18H 
I S MULTIPLICATION ) 

GO TO 132 

WRITE!N0UTPE,18 ) NM1,N 

F0RMAT!40H THE COMPOSITION OPERATOR BETWEEN STAGESI6,4H ANDI6,34H 



22 



aiE: DP4 



FORTRAN A1 NAVAL POSTSRADJATc SCHOOL 



IIS TO MAXIMIZE THE MINIMUM RETURN) 

GO TO 132 

14 WRITE(N0UTPE,19) NMl , N 

19 F0RMAT(40H THE COMPOSITION OPERATOR BETWEEN STA3£SI6,4H ANDI6,34H 
IIS TO MINIMIZE THE MAXIMUM RETURN) 

132 DO 136 I=l,ITOP 
134 BEST=1.0E+35*XM00E 
X = I-1 

XN=XLOW+X*0£LX 
NUFF=1 
CALL OLIMIT 
IF(MOOES) 602f602,603 

603 CALL S£ARCH(DLOWf DHIGHtDELD,XMOOE t DNBESTt BEST ) 

GO TO 400 

602 KTOP=( DHIGH-OLOW) /DELD+1.001 
223 DO 137 J=1,KT0P 
X=J- 1 

DN=0 LOW+X*DELD 

C JUMP MAY BE SET TO 1 IN TRANFM TO INDICATE AN INFEASIBLE XN-1 

C JUMP=0 

CALL TRANFM 
CALL STGRET 

C XM0DE=1 FOR MIN» -1 FOR MAX 

C I F( JUMP.E C.l ) RN=999999#XM0DE 

IF(NUFF) 400.401,401 

401 I F(N-1 ). 504,504 ,503 
504 QN=RN 

GO TO 143 

503 K=(YN-YLOW)/DELY+ 1.001 
IF(K) 212,212,209 
212 WRITE(N0UTPE,210 ) N,XN,DN,YN 

210 F0RMAT(///9H AT STAGE 15, 9H WITH XN=E15.8,8H AND 0N=£15.8, 12H XN-1 
1EQUALSE16.8,21H AND IS OUT OF LIMITS///) 

GO TO 2003 

209 IF(K-JTOP) 205,206,212 

206 RNM1=FNM1(JT0P) 

GO TO 207 

205 X=K-1 

X=YLOW+X*OELY 

RNM1 = FNMUK) + (FNMUK«-1)-FNM1(K)) * I YN-X)/DELY 

207 GO TO (21,22,23, 24, 248), NNFLAG 

248 WRITE( N0UTPE,249 ) N 

249 FORMAT(43H1COMPOSITIGN OPERATOR NOT DEFINED FOR STAGE, 15) 

GO TO 2003 

21 QN=RN+RNM1 
GO TO 143 

22 QN=RN=!'RNM1 
GO TO 143 

23 IF(RN-RNMl) 141,141,142 

141 QN=RN 

GO TO 143 

142 QN=PNM1 
GO TO 143 

24 IF(RN-RNMl) 145, 145,146 

145 QN=RNM1 
GO TO 143 

146 QN=RN 

143 I F{( GN-BEST)*XMODE) 144,137,137 

144 BEST=QN 
ONBEST=DN 

137 CONTINUE 

400 WRITE(NCALC) BEST,DNBEST 
136 CONTINUE 

REWIND NCALC 

WRIT E( NTAPE) XLOW ,XHI GH, DELX, ITD? , XM3DE 
DO 404 I = 1,IT0P 
404 READ (NCALC) DUM , 0 NOFX N ( I ) 

REWIND NCALC 

WRITE (NTAPE) (DNOFXN( II) ,I I=1,IT0P) 

DO 402 1=1, ITOP 

402 READ (NCALC) FN(I),DUM 
REWIND NCALC 
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YLOW=XLOW 
YHIGH=XHIGH 
DELY=D£LX 
JT0P = IT0P 
CONTINUE 

WRITE(NTAPE) (FN( 
END FILE MTAPE 
REWIND MTAPE 
READ(MTAPE) NMAX 
N=NMAX 
J=N-1 

IF(J) 407,407,201 
DO 522 1 = 1 , J 
READ (MTAPE) 
READ(MTAPE) 
READ(MTAPE) XLOW, 
READ (MTAPE) 
READ(MTAPE) ( FN ( I 
BACKSPACE MTAPE 
BACKSPACE MTAPE 
READ (INTAPE, 162 ) 
F0RMAT(2E15.8) 
NPAGE=NPAGE+1 
WRITE(N0UTPE,152) 
WRITE(NOUTP£, 2008 
FORMAT!/// • DATA 
IF(XMODE) 323,323 
WRITE(NOUTPE,325 ) 
F0RMAT(29H THE PR 
GO TO 327 
WRITE(N0UTPE,326 ) 
F0RMAT(29H THE PR 
WRIT S(NOUTPE, 323 ) 
FORMAT (42H XN IS 
ID XN = 1PE14.5) 

171 IF(XLOW-XNl) 175, 

176 WRITE(N0UTP5,177) 

177 FORMAT (47H0THE PR 
1 ANDE 16.8) 

GO TO 240 

175 IF(XN2-XHIGH) 181 
181 IF(XN1-XN2) 178,1 

178 I Xl= ( XNl-XLDw) /DE 
1X2= (XN2-XL0W )/DE 
6EST=1 .0E+35^'XM0D 
DO 182 J=IX1, 1X2 
IF(( FN(J)-BEST)*X 

1 33 8EST = FN( J ) 

JSAV E=J 
132 CONTINUE 

X= I A YF— 1 
X N=X LGW+X=«'DELX 
WRITE! NOUTPE, 184) 
164 F0RMAT(12H OPTIMA 
READ (MTAPE) (DNOF 
186 DN=DNOFXN(JSAVE ) 
CALL TRANFM 
204 WRITE(N0UTPE,188 ) 

188 F0RMAT(7H0 N1 

NN=0 

203 WRITE(N0UTPE,189) 

189 FORMAT! 17, 1P3E20. 
NN=NN+1 

N = N-1 

IF(N-l) 244,193,1 
244 GO TO 240 
193 DO 406 1=1,4 
406 BACKSPACE MTAPE 
READ(MTAPE) XLOW, 
READ(MTAPE) (DNOF 
L=( YN-XLOW)/DELX+ 
IF(L) 214,214,215 



FILE: 



123 

161 

197 



201 

522 

407 



162 



2008 

324 

325 

323 

326 

327 

328 
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II),II=1,IT0P) 



XHIG 

),I = 

XNl, 



NPA 
) XN 
CAR 
,324 
N 

06LE 

N 

OBLE 

XNl 

TO 



175, 

XLO 

OGRA 



,181 
73,1 
LX+1 
LX + 1 
E 

MODE 



XN, 
L XN 
XN(I 



8X2H 

N,X 

5) 



93 



XHIG 
XN( I 
1.00 



H,DELX, ITOP, XMOOE 
1 ,ITOP) 

XN2 



G E 

i, XN2 

D • , 2F20.3 /// ) 



M IS 


TO 


MINIMIZE 


AI6 


M IS 


TO 


MAXI MIZE 


AI6 


,XN2 








BE CHCOSEN OPTIM 


ALLY 



176 

W,XHIGH 

M ONLY HAS I NFORMATI 



,14H STAGE PROCESS) 

,14H STAGE PROCESS) 
BETWEEN XN=1PE14.5 

ON ON XN' BETWEENE16 



, 176 
76 

. 001 

.001 



) 1 83,182,132 



BEST 

= 1PE14.5,16H OPTIMAL RETURN= IPEl 4. 5) 
I ), I 1=1 , ITOP ) 



XNl 3X2HDN16X4HXN-1) 
N,DN ,YN 



H,DELX, ITOP, XMOOE 
I), I 1=1, ITOP ) 

1 



8H AN 
8 ,4H 
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214 L=1 

219 NP1=N+1 

WRITE(N0UTPE,210) NP1,XN,DN,YN 
NN=NN+7 

215 IF(L-ITOP) 218,216,217 

216 DN=0N0FXN(L) 

XN=YN 

GO TO 196 

217 L=ITOP 

GO TO 219 

218 XN=YN 
194 X=L-1 

X=XLOW+X*0£LX 

ON=DNOFXN{L J«-(ONOFXN( L+1 J -DNOFXN( L ) > * { XN-X ) / 0 ELX 
196 CALL TRANFM 

IF(NN-50> 203 ,203 ,220 

220 NPAGE=NPAGE+1 

WRITE( NOUTPE, 152 J NPAGE 
GO TO 204 

112 WRITE(NOUTPE, 113 ) MAX,ITCP 

113 F0RMAT(72H THIS PROGRAM LIMITS THE NUMBER OF DISCRETE STEPS OF THi 
ISTATE VARIBLE TOI 6,26H AND THIS PROBLEM REQJIRE5i6» 

2002 FORM AT Cl',' ERROR HALT') 

2003 WRITE( 6,2002) 

STOP 

2000 WRITE(6,2001) 

2001 FORMATCl',' END OF DATA FILE') 

STOP 

END 

CSEARCH DISCRETE FIBONACCI SEARCH SUBROUTINE 

SUBROUTINE SEARCH (AA, BB, DELYY, XXMOOE, YYBEST,BEST) 

DIMENSION F(150) 

OPTIMIZE WITH RESPECT TO Y BETWEEN AA AND BB IN STEPS OF DELYY 
STORE OPTIMUM Y IN YYBEST AND THE OPTIMUM VALUE OF THE OBJECTIVE 
FUNCTION IN BEST 
XXM0DE=-1 FOR MAXIMIZE 
XXM0DE=1 FOR MINIMIZE 

SUBROUTINE MUST BE INITIALIZED BY CALLING IT WITH DELYY=-1.0 AT 

LEAST ONCE BEFORE IT IS USED FOR A SEARCH 

DELY=OELYY 

A=AA 

B = B3 

XMODE=XXMODE 
IF(DELY) 100,100,101 

100 F ( 1) =1.0 
F (2) =2.0 

DO 1 1=3,150 
FU)=F( I-D+F ( 1-2 ) + l. 0 
I F<F ( I )-1.0E+35 ) 1,2,2 

2 I I = I 
RETURN 
1 CONTINUE 
I 1 = 1 50 
RETURN 

101 YHI=B 
YLO= A 

FNO= ( YHI-YL0)/DELY+1.0 
DO 5 N=1,II 
IF(FNO-F(N)) 6,6,5 

5 CONT INUE 
WRITE! 6,7) 

7 F0RMAT(38H1ERR0R, TOO MANY POINTS TO BE SEARCHED) 

2002 FORMATCl',' ERROR HALT') 

2003 WRITE! 6,2002) 

STOP 

6 IFIN-2) 42,41,40 
40 Y 1=F ! N-2 ) *DELY + YLO 

CALL FUNCTM!Y1,GY1) 

Y2=F !N-1) ^i'OELY+YLO 
CALL FUNCTN! Y2,GY2) 

16 N=N-1 
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IF(N-l) 102,102tl03 

103 I F( ( GY2-GY1)«XM0DE) 19,19,20 

19 YL0=Y1+DELY 
Y 1=Y 2 
GY1=GY2 

Y2=F (N-1)*0ELY+YLG 
IF(Y2-BI 104,104,105 
105 GY2=1-0E+35*XM00E 
GO TO 16 

104 CALL FUNCTN( Y2,GY2» 

GO TO 16 

20 YHI=Y2-DELY 
Y2=Y1 
GY2=GY1 

IF(N-2) 34,34,35 

34 FNM2=0.0 
GO TO 36 

35 FNM2=P(N-2) 

36 Y1=FNM2*0ELY+YL0 
CALL FUNCTNI Yl, GYl) 

GO TO 16 

102 IF( ( GY2“GY1)*XM00E) 110,111,111 

110 YYBEST=Y2 
BEST=GY2 
RETURN 

111 YYBEST=Y1 
BEST=GY1 
RETURN 

41 Y1=A 

CALL FUNCTNI Y1,GY1I 
Y2=6 

CALL FUNCTN(Y2 ,GY2) 

GO TO 102 

42 Y1=A 

CALL FUNCTN( Y1,GY1J 

GO TO 111 

END 

CFUNCTN 

SUBROUTINE FJNCTN(Y,GY) 

DIMENSION FN( 10001) ,FNM1 ( 10001 ), 0 NOFX N ( 1 0001 ) 

REAL * 8 TAB( 4) 

COMMON /CON/TAB 

COMMON XN,N,DN, YN , RN , SLO W , XHI GH, D E LX , OLOW , DH I GH, D ELD, NUMB ER 
common YLOW, YHIGH ,DELY ,NUFF,FN,JT0P,NNFLAG ,NMAX, XLAM 
COMMON A,B,C,0, E, SUMELT, KODE 
EQUIVALENCE ( FN ( 1 ) ,FNM1( 1 ) ,ONOFXN(D) 



0N=Y 

N0UTPE=6 

C JUMP MAY BE SET TO 1 IN TRANFM TO INDICATE AN INFEASIBLE XN-1 

C JUMP=0 

CALL TRANFM 
CALL STGRET 

C XM0DE=1 FOR MIN, -1 FOR MAX 

C IF( JUMP.EO.IJ RN = 959599*XM00£ 

IF(N-l) 146,146,401 
401 K = (YN-YLOW)/OELY+1.001 
IF(K) 212,212,209 
212 WRITE(N0UTPE,210) N,XN,DN,YN 

210 F0RMAT(///9H AT STAGE 15, 9H WITH XN=E15.3,8H AND DN=£15.8, 12H XN-i 
1EQUALSE16.8,21H AND IS OUT OF LIMITS///) 

2002 FORMATCl',* ERROR HALT') 

2003 WRITE(N0UTPE,2002) 

STOP 

209 IF(K-JTOP) 205,206,212 

206 RNM1=FNM1 ( JTDP) 

GO TO 207 

205 X=K-1 

X=YLOW+X>i=OELY 

RNM1 = FNM1 (K) + (FNM1(K'- 1 )-FNMl(K )) * ( YN-X l/DELY 

207 GO TO (21,22, 23,24) ,NNFLAG 
21 QN=RN+RNM1 

GO TO 137 
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22 QN=RN«RNM1 
GO TO 137 

23 IF(RN-RNMl) 141, 

141 ON=RN 

GO TO 137 

142 QN=RNM1 
GO TO 137 

24 IF(RN-RNM1J 145, 

145 QN=RNM1 
GO TO 137 

146 QN=RN 
137 GY=QN 

RETURN 

END 



A1 NAVAL POSTGRADUATE SCHOOL 
141 ,142 
145,146 



27 



FILE: STORE 



FORTRAN A1 NAVAL POST3RAOUATE SCHOOL 



SUBROUTINE STORE 

DIMENSION FN( 10001) fFNMl ( 1 0001 ), 0 NOFX N ( 1 0001 ) , A( 1 01 ) , BI 10 1 ) , C ( 101 J 
DIMENSION DdOl ) , E(lOl) 

REAL ♦ 8 TA8( 4) 

COMMON /CON/TAB 

COMMON XN,N,DN, YN ,RN, SLOW , XHIGH, 0 E LXi DLOW , DH I GH , D ELD , NUMB ER 
COMMON YLOW, YHIGH , DEL Y ,NU F F , FN t J T 0 P, N NFL AG , NM AX , X LAM 
COMMON A,B.CfO,E, SUMELT, K00EltK0DE2 
EQUIVALENCE ( FN( 1 ) , FNMII 1 ) ,DNOFXN( 1)) 

I NDATA=8 
N0UTP£=6 

READ (INDATA, 803) K0DE1,K00E2 
808 FORMAT (213) 

I F(KODEl.EQ.O) WR ITEI NOUTPE f 1 1 0) 

I F(KODEl.EQ.l ) WRITE! NOUTPE, 101) 

I F(K0DE1.EQ.2 ) WRITE! NOUTPE ,202) 

I F(K00E1.EQ.3 ) WRITE( NOUTPE, 303) 

110 FORMAT (30H THE OBJECTIVE IS COST ) 

101 FORMAT (30H THE OBJECTIVE IS SMA ) 

202 FORMAT (30H THE OBJECTIVE IS MS^T ) 

303 FORMAT (30H THE OBJECTIVE IS AVAILABILITY) 

I F(KODE2.EQ.O) WR ITEI NOUTPE, 900) 

I F(K0DE2.EQ.l ) WR IT E( NOU T PE , 90 1) 

I F(KODE2.EQ.2 ) WR I TE I NOU T PE , 902 ) 

I F(K0DE2. EQ.3 ) WR ITEI NOUTPE, 903) 

900 FORMAT (30H CONSTRAINT ON COST ) 

901 FORMAT (30H CONSTRAINT ON SMA ) 

902 FORMAT I30H CONSTRAINT ON MSRT ) 

903 FORMAT (30H CONSTRAINT ON AVAILA3ILITY ) 

DO 333 K=1,NMAX 

READ (INDATA, 999 ) A IK ) ,C I K ) , E I K) , B I K ) , 0( K ) 

999 FORMAT (5F10.4) 

WRITE! NOUTPE, 999) A ! K ) , B I K ) , C ! K) , 0 I K) ,E(K) 

333 CONTINUE 
SUME LT=0. 

DO 30 I=1,NMAX 

30 SUMELT = SUMELT<-E I I )=S‘A! I )»dI I) 

WRITE! 6,765)SUMELT 
755 FORM AT (F10.3) 

RETURN 

END 
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, STORE T 

SUBROUTINE STGRET 

DIMENSION FN{ 1000 IJ ,FNM1 ( 10001 ) , 0 NOFX N ( 1 0001 ) ,A( 1 01 » , 8 ( 10 1 J , C ( 101 ) 
DIMENSION DdOl ) , E( 101) 

REAL * 8 TAB(A) 

COMMON /CON/TAB 

COMMON XNtN.ON, YN ,RN, SLOW ,XHIGH, OELX, DLOW, DHI GH, DELD, NUMB ER 
COMMON YLOW.YHI GH ,DEL Y ,NU FF , FN , J T 0 P , 'I NFL AG , NM AX, X LAM 
COMMON A,B,C,0,E, SUME LT, KCDEl , KOD E2 
EQUI VALENCE ( FN ( 1 ) , FNMK 1 ) ,DNOFXN ( 1) ) 

I F(KODE1.EO.O )G0 TO 7 7 
IDN=INT(DN) 

A6=A (N )*B (N) 

I S=I DN 
TERM = 0. 

T EMP=0. 

I F( I S.LT.O) 30 TO 11 
TERM = EXP( -A6) 

TEMP=TERM 

IF(IS.EQ.O) 30 TO 11 

00 10 1=1, IS 
TEMP=TEMPxtAB/ I 
IF(TERM.GE.. 99999) GO TO 11 

C I F(TEMP.LE..00001 ) GO TO 11 

10 TERM = TERM +TEMP 

11 CDF=TERM 
20 CONTINUE 

1 F(KOOEl.NE.l )G0 TO <V0 
FOLLOWING CHANGED 120 133 

SMA= ( AE*( l.-TEMP ) + ( lON-AB ) «( 1 .-3 0 F ) )/ A3 
SMA=( l.-TEMP)+( ION-AB)*( l.-COF)/AB 
RN=E (N)*A6*SMA/SUM£LT 
RETURN 

40 TWUS = ( l.-COF) AB*AB-2.*AB*I0N+I 0N*( I ON+1) ) / ( 2. <=A ( N) ) 
X+TEMP*B(N)*(AB-I0N)/2. 

I F(K0DE1.EQ.3 )G0 TO 5 0 
RN=E (N )«TWUS/SUMELT 
RETURN 
50 CONTINUE 

AMSRT=TWUS/AB 
AMTB F = 1./A(N) 

AMTTR = D(N) 

RN=( AMTBF)/(AMTBF+AMTTR+AMSRT) 

RETURN 

77 RN=C(N)*DN 

C IF(K0DE2. EQ.l .AND.N.EQ.i . AND.YN.LT.XN )RN = 999999 

RETURN 
END 
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CDLIMIT 

SUBROUTINE DLIMIT 

DIMENSION FN( IDOOll.FNMU 1 030 1 J , D NOFX N ( 1 0001 J , A ( 1 01 i ♦ 3 ( 10 1 J t C ( 101 
DIMENSION D( 101 ) , El 101» 

REAL * 8 TA8(4i 
COMMON /CON/TAB 

COMMON XN,Nt ON. YN.RN. SLO W , XH I GH, 0 E LX. DLOW .OHI GH . 0 ELD , NUMB E R 

COMMON YLOW.YHIGH .DEL Y .NU FF , FN . J T OP .N NFL AG .NMAX. X LAM 

COMMON A.8.C.D.E , SUME LT. KCDcl .K0DE2 

EQUIVALENCE ( FN ( 1 J . FN M i( 1 ) .ONOFXN (1)1 

I F(KODE2.NE.O )G0 TO 3 7 

DL0W=0. 

0ELD = 1. 

DHIGH=AMIM(XN/C( N) .50.) 

RETURN 

37 AB=A(N)=i'B(N) 

DO 555 10=1.30 
IS=I 0-1 
TERM=EXP(-A3) 

TEMP=TERM 

IF(IS.EQ.O) GO TO 11 
DO 10 1=1 .IS 
TEMP=TEMP«AB/I 
IFITERM. GE.. 99999 ) GO TO 11 
C I FITEMP.LE..00001 ) GO TO 11 

10 TERM = TERM +TEMP 

11 CCF=TERM 

AO TWUS = ( l.-CDF) »( A8*A8-2.»AB=*‘IS+IS» ( I S*- 1 ) ) / ( 2. * A< N ) ) 
X+TEMP»B(N)=:'(A6-IS )/2. 

IF<<OOE2.EQ.2.ANO.E(N)*TWUS/SUMELT.LE .XN)GO TO 556 
C CHANGED FOLLOWING TO 555 FROM 556 
IFIK0DE2.EQ.2 )G0 TO 555 
I FIK0DE2.EQ.3 )G0 TO 554 
SSS= (l.-TEMP) +( IS-AB) *(1.-CJF) /A3 
RR=E (N)^AB*SSS/SUMELT 
Y=XN-RR 
OLOW=0. 

NUMBER IN THE NEXT LINE HAS A BIG EFFECT, IT SHOULD NOT. 
IF(N.EQ.1)DL0W=1000. 

I FlN.EC.l.ANO.y .GT.O. .AND. I0.LT.30 )G0 TO 555 
IFiN.EC.l.ANO.Y.LE, 0. )DLOW=IS 
DHIGH=30.+DL0W 

I FOLQW.E0.10 00. ) 0HIGH=10 00. 

DELD = 1, 

RETU RN 

AMSRT=TWUS/A6 
AMT8 F=1./A(N) 

AMTTR=D(N) 

A V=AMTBF/(AMTBF + AMTTR+AMSRT) 

I F(AV.GE.XN) GO TO 556 
CONT INUE 
I S=30, 

DLOW=IS 
DELD=1. 

DHIGH=50. 

RETURN 
END 

BLOCK DATA 
COMMON /CON/TAB 

REAL * 8 TAB14) /' SUM*,* MUL T ' . * M AXM I N* , ' MI NMAX* / 

END 



C THE 
136 

554 

555 

556 
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TRANFM 

SUBROUTINE TRANFM 

OIMEnIiON 

REAL *8 TAB(4J 
COMMON /CON/TAB 

COMMON XNtN.DN, YN f RNt SLOW , XHIGH, DELX, GLOW , OH i GH, D ELD , NUMB £R 
COMMON YLOW, YHI GH , DEL Y ,NU FF , FN , J T OP , N NFL AG ,N M AX , X LAM 
COMMON A,B,C»0,E, SUME LT , K CDE 1 , KOO E2 
EQUI VALENCE ( FN ( i ) , FNM1{ 1 ) .DNOFXN ( 1» » 

I F(K0DE2.NE.D »GO TO 12 

YN=XN-C(N )*ON 

RETURN 

12 IDN=INT(DN» 

AB=A{N»*B(N» 

I S=I DN 
TERM=0. 

TEMP=0. 

IF(I S.LT.O) GO TO 11 
TEPM = EXP<-AB» 

TEMP = T£RM 

IF(I S.EQ.O) GO TO 11 
DO 10 I=1»IS 
TEMP=TEMP*AB/I 
I F(TERM.GE..99999 ) GO TO 11 
; IFCTEMP. LE.. 00001 ) GO TO 11 

10 TERM = TERM +TEMP 

11 CDF=TERM 

40 TWUS = (1.-COF)«(AB*AB-2.#AB*IDN+ION*(IDN + 1» )/(2.*A(N) » 

X + TEMP*B(N »*(AB-I DN)/2 . 

IF(K0DE2*EQ.2» GO TO 50 
IF(KODE2.EQ.3 » GO TO 51 
SMA= ( 1 .-TEMP» +( IDN-AB )♦( l.-CDF)/AB 
YN=XN-E(N )*AB*SMA /SUMELT 



IF(YN.LT.O) YN=0.0 



I F(N .NE.l.AND.YN. LT.O . OR . N . NE . 1. A N D. Y N ,GT . 1 . J JUM P = 1 
IF(N.EQ.1.AND. YN. GT.O )JUMP=1 
IF(JUMP.EQ.1>YN = 0. 

RETURN 

51 AMSRT=TWUS/AB 
AMTBF=1./A(N) 

AMTTR = D(N» 

AV=( AMTBF>/( AMTBF+AMTTR+AMSRT) 

YN=XN/AV 
C AV IS AVAIL 

RETURN 

50 YN=XN-E(N »«TWUS/SUMELT 
I F(YN.L£.0. )YN=0. 

RETURN 

END 



C( 101 » 
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DIME NS ION A(20» , B (20) , C( 2 0 ) , D ( 20 ) , E ( 2 0 ) , M I NQ ( 20) 

DIMENSION T(1 0,30 ) ,F( 1200 ) ,DX( ID, 1200 ) ,TN( 10, 31 ) 

INTEGER BUDGET, BB ,BBP 1,DX 

REALMS TIME,CTIME 

INDATA=8 

N0UTPE=6 

READ( INDATA, 888 ) NMAX, BUDGET 
888 F0RMAT(2I5) 

READ (INDATA, 777) ( MI NQ( I ) , 1 = 1 ,NIM AX) 

777 FORMAT(20I5) 

DO 333 K=1,NMAX 

READ (INDATA, 999) A (K ) , C ( K ) , E ( K) , B ( K ) , D( K ) 

999 FORMAT (5F10. A) 

C WRITE(NCUTPE,999) A (K ) , B ( K ) , C ( K) , D ( K ) , E( K ) 

233 CONTINUE 

TIME=CTIME(1 ) 

SUMELT=0. 

DO 30 1=1, NMAX 

30 SUMELT = SUMELT+£(I )*A( I )#8( I) 



DO 10 1=1 , NMAX 
DX( I ,1 ) = 0 
DO 110 KP1 = 1,31 
K=KP1-1 
AB=A ( I )=<‘B( I ) 

TERM=EXP(-A6) 

TEMP=TERM 

IF(K.EC.O) G3 TO 11 
DO 14 JI = 1,K 
TEMP=TEMP*AB/ JI 

C I F(T ERM.GE.. 999999) GO TO 11 

C IF(TEMP, LE.. 0000001 ) GO TO 11 

14 TERM=TERM +TEMP 
11 CDF=TERM 

40 TWUS = ( l.-CDF) ^( AB*AB-2.*AB*K+K=S'( < + 1 ) ) / (2 . *A( I ) ) 
X + TEMP*B(I )*(AB-K) /2. 

110 TN(I ,KP1) = E( I j^t^TWUS/SUMELT 
DO 70 K=l,30 

70 T ( I , K) =TN( I ,K+1 )-TN( I ,K) 

10 CONTINUE 

C DO 456 IJK=1,10 

C 456 WRITE( 6,555)(TN(KI,IJK),KI = 1,NMAX) 

C 555 FORMAT (5F12.7 ) 



DO 80 1 = 1 ,NMAX 
LIM=30/MINQ( I ) 

DO 85 K = 1 ,LIM 
TMOD=0. 

MM=M INC( I ) 

DO 90 J=1,MM 

90 TMOD = TMOD+T( I , ( K- 1 ) INQ ( I ) + J ) 

T( I, K)=TMOD 
85 CONT INU5 

C( I) =MINQ(I )#C( I ) 

30 CONTINUE 

I BD=BUDG£T+1 
DO 100 BBP1=1,IBD 
BB=BBP1-1 
XMIN=0. 

I STAR=0 
TRY=0. 

DO 210 1=1, NMAX 
I F ( BB-C( I ) .LT.O) GO TO 210 
I F( BB-C( I ) .GT.O) GO TO 211 
TRY=AMIN1( TRY , T( I, 1 ) ) 

GO TO 212 
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211 TRY= T( I ,OX( I ,BBP1-C(I ) ) + l ) +F(BBP1-C( I ) ) 

212 IF ( TRY.GT.X'^INI GO TO 210 
XMIN = TRY 

I STAR = I 
210 CONT INUE 

F (BBPl )=XMIN 

19 DO 21 I=ltNMAX 

IF ( ISTAR.EQ.O.OR.BB-CdSTAR) .LE. 0 J DX(I,B8P1) = 0 
IF ( ISTAR.EQ.O.OR.BB-C(ISTARULE.O) GO TO 21 
OX (I TBBPi)=ox(i,aePi-cnsTAR) » 

21 CONTINUE 

IF(I STAR.EQ.OIGO TO 100 

IF(8 8Pl-C(ISTAR).LE.l > DX( 1ST AR, 3 6 +1 ) = 1 
IF(8BP1-C(ISTAR) . LE.l ) GO TO 100 

OXU STAR, BBPl ) = AM INOl 0X( ISTAR,BB? 1-C( ISTARI) +1,30» 

100 CONTINUE 
KLOW=I 8D-5 
KHIGH= IBD 

DO 101 KBD=KLOW,KHIGH 
T IME = T IME-CTI ME( 1 ) 

C WRITE(6,3456JTIME 
3456 FORMAT! • TIM£=* , FI 2.91 

WRIT El 6, 1000) (0X( I.KBO), I = 1,NMAX) 

WRIT El 6,1000) I MIN Cl I) ,1=1 ,NMAX) 

1000 F0RMAT110I5) 

AMSRT=0. 

00 88 1=1 ,NMAX 

88 AMSRT = AMSRT+TNl I , MINQ I I ) *0X1 I , K3 0 ) ♦ 1) 

KBDM1=KBD-1 

WRITE 16,2000) AMSRT,K3DM1 
2000 F0RMAT11X,'MSRT= *,F10.8,‘ BJOGET= ’,17) 

101 CONTINUE 
RETURN 
END 

FUNCTION CTIMEl I ) 

REALMS CTIME 
DATA IFL/0/ 

IFl IFL.NE.O ) GO TO 10 

IFL = 1 

CTIME = 0.000 

CALL SETIME 

RETURN 

10 CALL GETIMEIITIME ) 

CTIME = -ITIME*0. 00002600 

RETURN 

END 
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